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We have derived a semiclassical trace formula for the level 
density of the three-dimensional spheroidal cavity. To over- 
come the divergences occurring at bifurcations and in the 
spherical limit, the trace integrals over the action-angle vari- 
ables were performed using an improved stationary phase 
method. The resulting semiclassical level density oscilla- 
tions and shell-correction energies are in good agreement with 
quantum-mechanical results. We find that the bifurcations of 
some dominant short periodic orbits lead to an enhancement 
of the shell structure for "superdeformed" shapes related to 
those known from atomic nuclei. 
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Introduction — The periodic orbit theory (POT) 0-0] 
is a nice tool for studying the correspondence between 
classical and quantum mechanics and, in particular, the 
interplay of deterministic chaos and quantum-mechanical 
behavior. But also for systems with integrable or mixed 
classical dynamics, the POT leads to a deeper under- 
standing of the origin of shell structure in finite fermion 
systems from such different areas as nuclear P,p|-pl)[ , 
metallic cluster or mesoscopic semiconductor 

physics • Bifurcations of periodic orbits may have 

significant effects, e.g., in connection with the so-called 
"superdeformations" of atomic nuclei , and were 

recently shown to affect the quantum oscillations ob- 
served in the magneto-conductance of a mesoscopic de- 
vice 11. 

In the semiclassical trace formulae that connect the 
quantum-mechanical density of states with a sum over 
the periodic orbits of the classical system [Q-^, diver- 
gences arise at critical points where bifurcations of pe- 
riodic orbits occur or where symmetry breaking (or 
restoring) transitions take place. At these points the 
stationary-phase approximation, used in the semiclassi- 
cal evaluation of the trace integrals, breaks down. Vari- 
ous ways of avoiding these divergences have been studied 
[ 2|^ 16 ] , some of them employing uniform approximations 
[1^20|. Here we employ an improved stationary-phase 
method (ISPM) for the evaluation of the trace integrals 
in the phase-space representation, based on the studies 
in Rcfs. which we have derived for the elliptic bil- 

liard It yields a semiclassical level density that is 

regular at all bifurcation points of the short diameter or- 



bit (and its repetitions) and in the circular (disk) limit. 
Away from the critical points, our result reduces to the 
extended Gutzwiller trace formula (3|j|-|^ and is identi- 
cal to that of Berry and Tabor |^ for the leading-order 
families of periodic orbits. 

The main purpose of the present note is to report on 
the extension of our semiclassical approach to the three- 
dimensional (3D) spheroidal cavity which may be 
taken as a simple model for a large deformed nucleus 
or a (highly idealized) deformed metal cluster |l^,|l^ , and 
to specify the role of orbit bifurcations in the shell struc- 
ture responsible for the superdeformation. Although the 
spheroidal cavity is integrable (see, e.g., Ref. [^), it ex- 
hibits all the difficulties mentioned above (i.e., bifurca- 
tions and symmetry breaking) and therefore gives rise to 
an exemplary case study of a non-trivial 3D system. We 
apply the ISPM for the bifurcating orbits and succeed 
in reproducing the superdeformed shell structure by the 
POT, hereby observing a considerable enhancement of 
the shell-structure amplitude near the bifurcation points. 

Theory — The level density g{E) is obtained from the 
semiclassical Green function Q by taking the imaginary 
part of its trace in (1,0) action-angle variables 
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X exp I J [So. (I', I", i„) + (I" - I') • ©"] - *|A^a| . (1) 

Here {si} is the single-particle energy spectrum and H = 
H{T) is the classical Hamiltonian. The sum is taken over 
all classical trajectories a specified by the initial actions 
I' and final angles 0". (I', I", t„) = - J^" I -de is the 
action integral and ta the time for the motion along the 
trajectory a, and Ha is the Maslov index related to the 
caustic and the turning points | |2l] , p^ . In the spheroidal 
variables {u, v, ip}, the action I has the components 

lu — — J du \J (7i — sin^ u — (T2 / cos^ u , 

— lie 
Vt 

ly = — J dv \J cosh^ u — (Ti — (T2/ sinh^ u , 
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Hereby p = {2mEy/'^ is the particle's momentum and 
c = {b^ — a^y^^ half the distance between the foci; b and 
a are the semiaxes (with b > a) oi the spheroid with its 
volume fixed by a^b = and the axis ratio rj = b/a 
as deformation parameter; and ±Uc (or Vc) and vt are 
the caustic and turning points, respectively. In Eq. (||) 
we use the dimensionless "action" variables cti , a2 in 
which the torus of the classical motion is given by 

(To^ = < 172 < ^ 



if - 1 

,2 



172 < CTi < 



. . fT2 (77^ - 1) = a+. (3) 

In the ISPM, we expand the action Sa in up to 
second-order terms around its stationary points and keep 
the pre-exponential factor at zero order, taking the inte- 
grations over the torus within the finite limits given by 
Eq. (H). For the oscillating ("shell-correction") part of 
the level densit y S g(E) = g{E) ~g{E), where g{E) is its 
smooth part [0j2J| , we obtain 



Sg{E) 



J- Re ^ {E) exp (^ikLp 



where k = p/fi is the wave number and Eq = ff /2mE? 
our energy unit. The amplitudes Ap will be specified 
below. The sum over /? includes all 2-parameter fami- 
lies of three-dimensional (3D) periodic orbits and ellip- 
tic and hyperbolic 2D orbits lying in a plane contain- 
ing the symmetry axis (all with degeneracy parameter 
/C=2), the 1-parameter families of (2D) equatorial orbits 
lying in the central plane perpendicular to the symme- 
try axis (with /C=l), and the (ID) isolated long diam- 
eter (with /C=0). In Eq. (^, Lp is the length of the 
orbit (3 at the stationary point (o'^,(T2) which for the 
3D orbits lies inside the physical region of the torus (||), 
and is analytically continued outside this region. The 
(T2 = boundary of (^) is occupied by the 2D orbits 
with IC—2. The stationary points are determined by the 
roots of the periodicity conditions uju/^v = nu/n-u and 
uju/w^p = riu/n^; hereby — 811/81^ are the frequen- 
cies and rifj are co-prime integers which specify the peri- 
odic orbits f3 = M{ny, n^, Uu), where M is the repetition 
number. The factor wj = exp(— 7^L^/4i?^) in Eq. (^) 
is the result of a convolution of the level density with 
a Gaussian function over a range 7 in the variable kR. 
This ensures the convergence of the POT sum (^) by sup- 
pressing the longer orbits which are not relevant for the 
coarse-grained gross-shell structure . 

For Strutinsky's shell-correction energy 5U |^,0,|4| , we 
obtain (with time reversal symmetry and a spin factor 2) 
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The Fermi energies Ep (and with it kp) and Ep are 
determined by the particle number conservation N = 

2 J^''g{E)dE = 2 J^''g{E)dE. Due to the factor 
the PO sum in ^ may converge faster for the shortest 
orbits than the level density (^ for small 7. Any en- 
hancement of the amplitudes Ap of the most degenerate 
short periodic orbits - e.g., due to bifurcations or to sym- 
metry restoring, as discussed below - therefore leads to 
an enhancement of the shell structure and hence to an 
increased stability of the system. 

We present here only the amplitudes of the leading con- 
tributions to (|^) and (||). For further details (including, 
e.g., explicit expressions for the Maslov indices /Xq), we 
refer to a forthcoming, more extensive publication [ p2[ . 

For the amplitudes Ap of the most degenerate {IC—2) 
families of periodic 3D and 2D orbits, we obtain 
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(4) The quantity 



K'jp k'^P is related to the main cur- 
vatures k'p^'^ of the energy surface E — H{ai, (72) in the 
"action" plane {(Ti,a2), given by 



In Eq. 



9(72 Q^2 g^2 



(n-1,2) (7) 



the arguments of the two-dimensional error 



function erf (a;, y) — 2 dz e 
turning points in the action plane 



TT are given by the 



= ^-mMn,K^;^lh [at - ; (n = l,2) (8) 

see Eq. (^ for the boundaries ct^. All quantities in (||) 
can be expressed analytically in terms of elliptic integrals. 
For the 3D orbits, our result (|) is in agreement with 
that obtained by exact Poisson summation over the EBK 
spectrum (cf. Refs. [||@]). 

For the contribution of the K, — 1 families of equatorial 
orbits to (0), we obtain the amplitudes 



A 



(K=l) 
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where 0^ — irn^/ny, Fp is their stability factor 
al = = cos^ 'Pp/iv'^ - 1), and 



X Q — kc 
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2:3=0. 



(10) 



The contribution of the isolated long diameter orbit, 
which may be expressed in terms of incomplete Airy in- 
tegrals [21 2^, is not important for deformations of the 
order -q ~ 1.2 — 2. 
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FIG. 1. Moduli of amplitudes |j4^| versus 77 for the equato- 
rial "star" orbit (5,2) {K. = 1) and the 3D orbit (5,2,1) {K. = 2) 
bifurcating from it at r;=1.618... Solid lines: using the ISPM 
according to Eqs. ^ and (^, respectively; dash-dotted lines: 
using the standard stationary-phase approach. 

Discussion of results — In Fig. 1 we show \Af}\ versus 
deformation 77 (at kR = 40) for a pair of orbits involved 
in a typical bifurcation scenario. At the critical point 
77 = 1.618... the equatorial "star" orbit (5,2) undergoes 
a bifurcation at which the 3D orbit (5,2,1) is born; the 
latter does not exist below 77 = 1.618... 

In the standard stationary-phase approach (SSPM; 
dash-dotted lines), the amplitude of the (5,2) orbit di- 
verges at 77 = 1.618..., whereas that of the bifurcated 
orbit (5,2,1) is finite but discontinous. As seen in Fig. 1, 
the ISPM (solid lines) leads to a finite amplitude 
for the (5,2) orbit. This is because the factor Fp in the 
denominator of (^, which goes to zero at the bifurcation, 
is cancelled by the same factor in the numerator of 
( |l0| ) via the third error function in (^. A similar result 
was found for the short diameter orbit 2(2,1) in the el- 
liptic bilHard Furthermore, the ISPM softens the 
discontinuity for the (5,2,1) orbit, leading to a maximum 
amplitude slightly above the critical deformation. 

The relative enhancement of these amplitudes Ap near 
the bifurcation point can also be understood qualitatively 
from the following argument. At the bifurcation of the 
equatorial (5,2) orbit, its degeneracy parameter /C = 1 
locally increases to 2, because it is there degenerate with 
the orbit family (5,2,1) that has K. — 2 at all deforma- 
tions 77 > 1.618... This is similar to a symmetry restoring 
transition. An increase of the symmetry parameter /C by 
one unit leads to one more exact integration compared to 
the SSPM, and thus the amplitud es (P) and (^) acquire 
an enhancement factor ^/kLp (x yJpR/Ti (cf. Refs. |§,0|). 

A similar enhancement of the double triangle 2(3,1) 
and the 3D orbit (6,2,1) is found near their bifurcation 
point 77 = VS = 1.732... However, the curvature k'"^^ (|^) 
for orbits like M(3t,t, 1) {t = 2,3,...) is identically zero 



and hence the SSPM is divergent for all deformations 77 > 
1, in contrast to the situation with orbits like (5,2,1) with 
finite K^jP ■ Here we have to take into account the next 
nonzero 3rd-order terms in the expansion of S^, although 
the (3t,t, 1) ISPM amplitude (||) is finite and continuous 
everywhere. The amplitude can then be expressed in 
terms of incomplete Airy and Gairy integrals with finite 
limits [|2^. For the equatorial orbits t(3,l), like for the 
double triangles 2(3,1), one has a zero curvature k'^^ 

only at the bifurcation point 77 = \/3- Here Ffj j 0, 
and a similar mechanism of cancellation of singularities 
for other orbits takes place through Eqs. ([^|l0|). But the 
relative enhancement of the ISPM amplitudes (||j9|) of 
such orbits at the bifurcations is of order fcL^ because 
of a change of the degeneracy parameter /C by two units 
(see for details). In this sense we avoid here a double 
singularity related to a double restoring of symmetry. 

In Figs. 2 and 3, we present semiclassical level densi- 
ties Sg{E) (^) versus kR and shell-correction energies SU 
versus A^^/"^ for various critical deformations (heavy 
dotted lines), and compare them to the corresponding 
quantum- mechanical results (thin solid lines). We ob- 
serve a very good agreement of the gross-shell structure 
at all deformations. The most significant contributions 
to these results near the critical deformations are coming 
from bifurcating orbits with lengths smaller than about 
lOR, in line with the convergence arguments for the POT 
sums (Q) and (^) mentioned above. For the bifurcation 
at 77 = 1.618..., the orbits (5,2,1) and (5,2) give contribu- 
tions comparable with other 2D orbits. For 77 = v^, the 
bifurcating orbits (6,2,1) and (6,2) are also important. 
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FIG. 2. Level density SgiE) (|) (unit E^'^) versus kR for 
different critical deformations 77. The Gaussian averaging pa- 
rameter is 7 = 0.3. Thin solid lines: quantum-mechanical re- 
sults; heavy dotted lines: semiclassical results using the ISPM. 



3 



200.0 I ' ^ ' ^ ' 

100.0 - Tl=1.2 J 

0.0 -—v, — ^-■^^''-vA./%A/\/AK^n/\^AAi 

-100.0 - 
-200.0 



-SUo, 



-8U,s 



100.0 

0.0 
-100.0 
-200.0 



ri=1.618... 



-5U„, 



•■■8U,s 



100.0 - ri=1.732... 

0.0 ----^■-~-'=*i;'-^'W-^^/VV>VVYV 
-100.0 - 



-5U, 



SUisPM ^ 



-200.0 I 

100.0 - 11=2 

0.0 —'-»-"-*"-V^i.f'VV- 

-100.0 - 

-200.0 

3.0 



-SUo, 



■■■SUis 



FIG. 3. Shell-correction energy 5U (g) (unit Eo) versus 
cube root of particle number A''^^'' (same notation and same 
deformations as in Fig. 2). 

The role of the bifurcating orbits increases for larger 
deformations and is dominating at the superdeformation 
r] = 2. For this deformation, the most important orbits 
in the present spheroidal cavity model are the 3D orbits 
(5,2,1), (6,2,1), (7,2,1), and (8,2,1). 

These results are in agreement with both heights and 
positions of the peaks in the length spectra obtained in 
Ref. from the Fourier transforms of the quantum level 
densities g{kR) at the same deformations. 

Summary and conclusions — We have obtained an 
analytical trace formula for the 3D spheroidal cavity 
model, which is continuous through all critical deforma- 
tions where bifurcations of periodic orbits occur. We find 
an enhancement of the amplitudes \Afj\ at deformations 
r] ~ 1.6 — 2.0 due to bifurcations of 3D orbits from the 
shortest 2D orbits. We believe that this is an important 
mechanism which contributes to the stability of superde- 
formed systems. Our semiclassical analysis may therefore 
lead to a deeper understanding of shell structure effects in 
superdeformed fermionic systems - not only in nuclei or 
metal clusters but also, e.g., in deformed semiconductor 
quantum dots whose conductance and magnetic suscep- 
tibilities are significantly modified by shell effects. 
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